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Abstract 

We extend Wegner's exact solution for the 2D density of states at the lowest 
Landau level with a short-range disorder to the cases of a double-layer system 
and a superlattice. For the double-layer system, an analytical expression for 
the density of states, illustrating the interplay between the tunnel splitting of 
Landau levels and the disorder-induced broadening, is obtained. For the su- 
perlattice, we derive an integral equation, the eigenvalue of which determines 
the exact density of states. By solving this equation numerically, we trace the 

disappearance of the miniband with increasing disorder. 
PACS numbers: 73.20.Dx, 73.40.Hm 
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The shape of the Landau levels (LL) in a 2D system in the presence of a disorder was the 
subject of intensive study during the last two decadeslHll The complexity of the problem 
arises from the fact that in the absence of the disorder the energy spectrum is discrete. 
As a result, the self-energy of an electron appears to be real in any finite order of the 
perturbation theory. Therefore, obtaining a finite width of the LL requires summation of 
the entire diagram expansion. It was demonstrated^'^ that such a summation is possible 
when the number of the LL is large. The simplifications, arising in this limit, are different 
in the case of a short-range and a smooth disorder. In the former case only a subsequence 
of diagrams without self-intersections contributes to the self-energy, or, in other words, the 
self-consistent Born approximation^! becomes asymptotically exact .0 The shape of the LL 
in this case is close to semielliptical. For a smooth disorder, with correlation radius larger 
than the magnetic length, all diagrams are of the same order of magnitude, but in this case 
magnetic phases, caused by self-intersections of impurity lines, become small. The origin of 
these phases lies in an uncertainty in the position of the center of the Larmour orbit. Having 
the phases dropped, the entire perturbation series can be summed up with the help of the 
Ward identity, resulting in the Gaussian shape of the LL.0 

For low LL numbers and short-range disorder, the magnetic phases in diagrams are of 
the order of unity. A small parameter appears in the problem only if the energy e (measured 
from the lowest LL) is much larger than the LL width T, making possible a calculation of 
the density of states (DOS) in the tails of LL. Such calculations were carried out in the 
framework of the instanton approachiffl'0ll!@ll! and the tails were shown to be Gaussian. 
In the domain e ~ V the problem has no small parameter and no simplifications are possible. 
However, for the lowest LL, the exact DOS was found by Wegneri for an arbitrary ratio 
e/Y. Wegner has shown that the diagrammatic expansion of the disorder-averaged Green 
function, G(e), can be mapped onto that of the zero-dimensional complex <£> 4 -model with 



the partition function given by a simple integral 



Zq~\s, Y) = J dip* dip exp 



Y 2 

ieip*ip - — (p*p) 2 
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The crucial observation made by Wegner was that the number of diagrams for the disordered 
system, which are mapped onto a single graph of the y9 4 -model, equals (up to an overall 
factor) the inverse value of the diagram itself. The electron Green function is then given by 

°~ 

where / is the magnetic length. Wegner has proved that coefficients in front of r n in each 
side of this equation coincide. Having a closed expression for G(e), Wegner obtained the 
following formula for the DOS in the lowest LL 

The magnetic field dependence of the width r is T oc \f~B. More precisely, for the correlator 
of the random potential V(r) of the form (V(r)V(r')) = 7<5(r — r'), one has T = (7 /2-kI 2 ) 1 ^ 2 . 
An alternative derivation of Wegner's result was given by Brezin, Gross and Itzyksonl in 
the framework of functional-integral approach. 

Consider now a system consisting of two parallel two-dimensional layers. In the absence 
of a disorder and magnetic field, a tunnel coupling between the layers would cause a splitting 
of size quantization levels by an amount of 2t, t being the tunnel integral. In a perpendicular 
magnetic field, the spectrum of the system represents two staircases of LL shifted in energy 
by 2t. Assume that the field is strong, so that the cyclotron energy is much larger than t. If 
a disorder is present in the layers, the resulting shape of two adjacent LL's would depend on 
the ratio T/t. If this ratio is large, then the tunnel coupling does not play any role, so that 
the DOS is twice the DOS in an individual layer, which is given by Eq. (|3|). In the opposite 
case, f > T, the peaks in the DOS, corresponding to the symmetric and antisymmetric 
combinations of size-quantization wave functions, are broadened independently. The centers 
of the peaks are distanced by 2t and their shape is described by Eq. (|3]) with the widthS 
T/y/2. The factor l/y/2 appears because the effective random potential for symmetric 
(antisymmetric) state is [Vi(r) ± V2(r)]/2, where Vi(r) and V^(r) are the random potentials 
in the layers. If (Vi(r)V2(r')) = 0, the correlator for each effective potential appears to be 



twice as small as that for an individual layer. 

In the case r ~ t, calculation of the DOS in a double-layer system seems to pose even 
harder problem than for a single layer, since here the DOS represents a two-parametric 
function, gdi(e/T,t/T), with both arguments of the order of unity Nevertheless, as we 
demonstrate below, for the lowest LL the exact DOS can be obtained in a closed form by 
generalizing Wegner's approach. Moreover, such a generalization can be carried out for an 
arbitrary number of layers, and, in particular, we consider the case when the number of layers 
is infinite (superlattice). In the absence of a disorder, each LL in a superlattice gives rise to 
a miniband of a width 4t. Gradual switching on a disorder first smears out the singularities 
in DOS at the edges of the miniband and then, as T exceeds t, transforms the DOS into 
a single peak, corresponding to an individual layer. We derive an integral equation, the 
eigenvalue of which determines the DOS in a superlattice, and trace this transformation by 
solving it numerically. 

Consider first the double-layer system. The Hamiltonian has the form 



+ ip*Viipi 



(4) 



where A = (— By/2, Bx/2) is the vector-potential in the symmetric gauge measured from 
the origin in both layers. It will be convenient to include the last term in (Q) into the 
definition of the free Hamiltonian. Then, after projecting onto the lowest LL, the free Green 
function represents a 2 x 2 matrix 



G°(v,v') 



_Q_ 

2tt/- 



exp 



J\2 



4/ 2 



with 
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The perturbation expansion of the Green function, averaged over random potentials V\ 
and V2, has the same diagrammatic representation as for a single layer. Some of the first 
diagrams are shown in Fig. 1(a) and (b). The solid lines correspond to G° and the dashed 



lines correspond to the correlator of the random potential. In contrast to the single-layer 
case, solid lines carry indices, reflecting the fact that electron can tunnel from one layer to 
another between two successive scattering acts. Since the scattering retains an electron in 
the same layer, the indices at the ends of each dashed line are the same. Let us introduce 
projecting operators fj as 



(l 
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Then the expressions corresponding to diagrams (a) and (b) can be written as 

#1) = F^T 2 Y,mQf l Q), 

i 

Q{2) = F (2) T 4 Y J {Qr i QT j QT i QT j Q), (8) 

ij 

where F^ and are spatial integrals. Similarly, in any nth order diagram the spatial 
integrals are separated out as factors in front of products of matrices Q and fj, which are 
responsible for the energy dependence. Important is that coefficients F^ are exactly the 
same as those for a single-layer. 

The mapping is carried out following Wegner's prescription: one identifies pairs of points 
in a diagram connected by dashed lines, and one gets a graph with four lines entering each 
vertex [see Fig. 1(c) and (d)]. In doing so, one obtains, in general, a set of diagrams yielding 
the same graph. It is clear, however, that since assigning indices to the lines does not alter 
in any way the topology of diagrams or graphs, the number of diagrams in a set is the same 
for both single- and double-layer cases. Moreover, one observes that the contractions of 
matrices fj precisely follow the identification of points described above [as it can be seen, 
e.g., in Fig. 1(b) and (d)], so that all the diagrams in such a set are equal. The fundamental 
relation, established by Wegner, is that for each diagram in the set one has F^ = 1/Afs, 
where Af is the number of diagrams in the set and 1/s is the symmetry factor of the graph (s 
is the number of permutations leaving graph invariant). The latter factor is also unchanged 
by assigning indices to the graph. For example, the graph (d) is invariant under permutation 
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of upper and lower lines so its symmetry factor is 1/2 in both cases. Thus, the contribution 
of the set, being proportional to J\fF^ n \ is TV-independent, and the problem again reduces 
to the zero- dimensional field theory. The remaining question is whether matrix products of 
type @ can be generated in the perturbation expansion of some generalized <£> 4 -model. Our 
main observation is that the model with the partition function 



Zf = J d$*d$exp 



2 



(9) 



accomplishes this task. Here Q and fj are matrices defined by @, and $ is a two- component 
complex field: $ = (^1,^2). Indeed, the nth order term in the expansion of exponent 
(§) in terms of T 2 represents a product of 2n matrices fj (with all pairwise contractions) 
separated by 2n products of the form $$*. Then the gaussian integral over $ inserts the 
"Green function" Q = in place of each pair of fields $ and $*, with all possible 

contractions between them yielding all the nth order graphs with appropriate symmetry 
factors. 

Having the mapping established, the DOS in the double-layer system can be calculated 
directly from (H). It is also instructive to rewrite Zq in a different form. First, we decouple 
the quartic term in the exponent of (H), 

r 2 



iS= £ 

i=l,2 L 



-itfalw + tp&i), (10) 



with the help of gaussian integral over a pair of auxiliary variables. Performing the remaining 
integral over (pi we then obtain 

zf = Mrr , , ,^ A v „ ex P (-M - . (id 



71T 2 J-00J-00 (e + Ai)(e + A 2 ) - t 2 ^ \ T 2 T 2 
From the form (|TTj) , the both limiting cases of large and small t are evident. For small t, 
the partition function factorizes, = (^o^) , yielding twice the DOS (0). For t 3> T 
the characteristic values of Ai, A2 in (0), being of the order of T, are much smaller than 
t. This allows to neglect the product A1A2 in the denominator; Zq is not small only if 
(e — t) ~ T or (e + t) ~ Y. In both cases one should introduce new variables /ii = Ai + A 2 



and H2 — Ai — A2. Then the integration over /i 2 would contribute a factor \/27rr, and the 
integral over fi x would reproduce Wegner's result with the width T/^/2, as discussed above. 
The evolution of DOS between two limits, calculated from (|2|) with Z$ 2 \ is shown in Fig. 2. 

Let us now turn to a superlattice. The partition function (|9|) can be straightforwardly 
generalized to a multilayer system, and for n layers with nearest-neighbor tunneling it takes 
the form 

n n—1 p2 n 

- H^2{<Pj<Pj+i + ¥>j+i¥>j) - -r XX^>i) 2 ■ (12) 
j=i j=i 4 3=1 J 

We are interested in the asymptotic behavior of Zq (e, T) as n — > 00. For this purpose we 
employ a method similar to the transfer-matrix method in the theory of ID spin chains. 
Note that the expression ( fL2D for Zq can be rewritten as 

r 2 



/n 
i=l 



zt\e,T) 



dip* dip exp 



i £( p*(p - ^-(ip*ip) 2 



In(<P*,<f), 



(13) 



where I± = 1 and the functions I n (}P*-> V 9 ) satisfy the following recurrence relation 

r 2 



I n +i(</, V?) = T e pl n = J dip* difx exp 



- ^(V^i) 2 - it(<p*<pi + (flip) 



WnVi). (14) 



Consider now the eigenvalues, X^ k \e, V), and eigenfunctions, ^^(ip* , ip) , of the operator 
T £j r: T £) rf2( fc ^ = X^fl^. Assume that A 1 - ** has the maximal absolute value. Then in the 
limit n — > 00, Zq™^ will behave as (A 1 - ))" - . Hence, the DOS per layer in a superlattice can be 
expressed through A^(e, T) in the following way 



9*i(e, T) 



- 1 1 1 1 lim ! ' 



1 T 91nA(°)(e,r) 
-lm- 



2tt 2 Z 2 n-+oo n Q £ 2tt 2 1 2 " ' de 
Thus, we have reduced the calculation of DOS to the solution of an integral equation 



(15) 



X {0 \e,T)Q { ^{ v *, v ) = J d<p* x dip x exv 



iecplcp! - — (y?>i) 2 - it((p*(pt + (p* x (p) 



rf&fart). (16) 



Consider first the case of a weak disorder, r — > 0. One can check that eigenfunctions of T e o 
in this case have the form 



(17) 



where R and a are, respectively, the absolute value and the phase of ip, and L™(x) is the 
Laguerre polynomial. The corresponding eigenvalues, \( p ' m \s, 0), are equal to 

\ 2p+m+l 



XM(s,0) = -( — 
t \ie 



2it 



V4t 2 - e 2 



(18) 



where for \e\ > 2t the square root is defined as i~ 1 sgn(e)^/e 2 — At 2 . Outside the interval 
|e| < 2t, the phases of eigenvalues (|I8| ) have no energy dependence, supporting the obvious 
observation that the DOS is zero outside the miniband. Within the miniband, all X (j> ' m \e, 0) 
have the same absolute value. This is a manifestation of the fact that for a large but finite 
number of layers the DOS in the absence of disorder represents a set of delta-peaks. However, 
with arbitrary weak disorder present, only the eigenvalue \^°'°\e, 0) will survive in the limit 
n — > oo, yielding the familiar result 

1 



1 dlnA^M) 
ga{£, 0) = ^^Im- 



(19) 



2tt 2 / 2 "" de 27r 2 /V4t 2 - e 2 ' 

Assume now that the disorder is finite but f < i. It is convenient to formally rewrite 



Eq. ( |l6l) in the following form 



(o) 



dE exp 



(E 



(20) 



y/TTT J-c 

For small T, only the energies E close to e contribute to the integral (p0|). This suggests 
to start the iteration procedure by substituting, as a zero approximation, the m = p = 
eigenfunction of Te,q, ^eo 



exp 



E + V4t 2 - E 2 ) R 2 /2 , into the right-hand side. 



This generates the first approximation for the function f^°p 



n 



(0) 



0FrA(°)( £ ,r) J- 



dEexp 



(E-e) 
T 2 



\(°>°\E,o)nf$ 



(21) 



Substituting this function back into (|20|), we obtain 



7rr 2 A(°)(e, T) J- 



dE / dE'exp 



(E 



(E> - e) 



r 2 



\^(E',o)f E ^;°l (22) 



Note now, that ^l^o as a function of E' changes significantly on the scale E' ~ t. On the 
other hand, exponential factors in (|22"D enforce the difference between E and E 1 be of the 
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order of T. This allows to replace fl^'o ^e'o under the integral. Then we immediately 



observe that the right-hand side takes the form X^(e, r)f2fp with 



A ( %,r) 



2tt r 00 dE 



exp 



(E-e) 

r 2 



(23) 



V^T J-oo V4t 2 - E 2 - iE 

In other words, for small Y the function Vtfl satisfies Eq. (|20|), yielding the eigenvalue (p3|). 

In principle, to assess the region of large T one should keep iterating Eq. (p0|) . However, 
as we have established numerically, the function fr ** is already a very good approximation 
for and X^ is a very good approximation for A^ not only for small, but for arbitrary 
ratio Y/t. The reason for this is the following. As Y/t increases, one should reproduce 
Wegner's result, which corresponds to t = and, hence, f^°) = const in Eq. fllTf ). On the 
other hand, it is easy to see that turns to constant as t — >• 0, and that in this limit X^ 
turns to Thus, Eq. ([23|) is exact in both limits. The numerical results for the DOS, 

obtained by substituting A^ -* into (15), are shown in Fig. 3. We see that the miniband is 
completely destroyed as Y/t exceeds 1.4. 

Note in conclusion, that a decade ago there was a significant interest in the study of 
transport phenomena in multilayer systems in a strong magnetic field (see, e.g., Refs. ^1] 
and |22p. In the recent publicationsi^E this interest was renewed. The question of interest is 
how a transition from a purely 2D to three-dimensional behavior of the conductivity occurs 
with increasing t. As was shown in Ref. the structure of electronic states in a multilayer 
system can be efficiently tuned by tilting magnetic field. In this case the role of the parallel 
component of B reduces to the suppression of the interlayer tunneling. 
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FIGURES 

FIG. 1. First (a) and second (b) order diagrams for the Green function mapped on graphs (c) 
and (d), respectively. 

FIG. 2. The DOS per layer for a double-layer system in units of g\ = (2tt1 2 )^ 1 T^ 1 for values 
of t/r=0.0 (highest curve), 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, and 3.0. 

FIG. 3. The DOS per layer for a superlattice in units of 52 = (27rZ 2 ) _1 (2i) _1 for values of 
r/2i=0.1 (highest curve), 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0. 
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